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Cepheids in Open Clusters: An 8-D All-sky Census 
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ABSTRACT 

Cepheids in open clusters (cluster Cepheids: CCs) are of great importance as zero- 
point calibrators of the Galactic Cepheid period- luminosity relationship (PLR). 

We perform an 8-dimensional all-sky census that aims to identify new bona- 
fide CCs and provide a ranking of membership confidence for known CC candidates 
through membership probabilities. The probabilities are computed for combinations of 
known Galactic open clusters and classical Cepheid candidates, based on spatial, kine- 
matic, and population-specific membership constraints. Data employed in this analysis 
are taken largely from published literature and supplemented by a year-round observ- 
ing program on both hemispheres dedicated to determining systemic radial velocities 
of Cepheids. 

In total, we find 13 bona-fide CCs, 3 of which are identified for the first time, 
including an overtone-Cepheid member in NGC 129. Inconclusive cases are discussed 
in detail, some of which have been previously mentioned in the literature. Our results 
are inconsistent with membership for 7 candidates that have been studied previously. 
We employ our bona-fide CC sample to revisit the Galactic PLR and obtain results 
consistent with most other calibrations, being limited by cluster uncertainties. 

In the near future, Gaia will enable our study to be carried out in much greater 
detail and accuracy, thanks to data homogeneity and greater levels of completeness. 

Key words: methods: data analysis, catalogs, astronomical data bases: miscella- 
neous, stars: variables: Cepheids, open clusters and associations: general, distance 
scale 



1 INTRODUCTION 

The search for Cepheids in Galactic open clusters (CCs) 
has been a topic of interest in astronomy for the past 60 
years, owing largely to their importance as calibrators of 
the Cepheid period- luminosity relation (PLR), discovered a 
century ago among 25 peri odic variable stars in the SMC by 
iLeavitt fc Pickerinel (|l912D . 

The proportionality between the logarithm of Cepheid 
pulsation periods and their absolute magnitudes, i.e., their 
luminosities, gives access to distance determinations and has 
established period- luminosity relationships as cornerstones 
of the astronomical distance scale (e.g. iFreedman et al.l 
I2OOII : ISandaee et aD I2OO6I). For re views on Cepheids as 
distan ce indicators, cf. FeastI () 19991 ): rSandage"fc Tammann 
l|2006h . for instance. 

The existence of the Cepheid PLR is most ob- 
vious among Cepheids in the Magel l anic Clouds (e.g. 
lUdalski et al.l 19991 : [Sos zvnski ct aL'2008.'2010l). due to com- 
mon distances (small dispersion), large abundance (thou- 
sands), and relative proximity (detectability). However, 



E-mail: r ichard . ander sonOunige . ch 



knowledge of the zero-point (s) of such relations is also re- 
quired; in this case, the distances to the Magellanic Clouds. 
For such zero-point calibrations, PLR-independent distance 
estimates are requi r ed, e. g. from trigonometric parallaxes 
iFeast fc Catchpolj Il997l : [Benedict et all I2007I'). Baade- 
Wess elink-type methods ijCieren et al.l 19971 : Storm et al.l 
l201lh . or o bjects l ocated at comparable dist ance, e.g. wa- 
ter m asers (|Macri et al. 2006 ') or open clusters (|Turner et al.l 
|2010D . 

For open clusters, distances can be determined via zero- 
age Main Sequence or isochrone fitting. If membership can 
be assumed at high confidence, the cluster provides the in- 
dependent estimation of the Cepheid's distance. Confidence 
in cluster membership is thus critical for such calibrations. 

Since the first discovery of CCs by llrwini (|l955l . iden- 
tified S Nor in NGC 6087 and U Sgr in M 25) and iFeastI 
(| 19571 . established membership via radial ve locities), many 
researchers have contribut ed to this field, e.g . Ivan den Berghl 
(1957); Efre moy (1964 ) : iTsareyskv et al.i lll966l) : iTurnp 



(il98d): iTurner et al.l ll993h: iBaumgardt et al 



Hovle et all (l2003l l: lAn et al.l (|2007l ): iMaiaess erall 



200C) 



200s 



Turneil (|2010h 7 Nevertheless, relatively few bona-fide CCs 



(< 30) have thus far been discovered. 



2 R.I. Anderson et al. 



We therefore carry out an all-sky census of classical 
Cepheids in Galactic Open Clusters that aims to increase 
the number of bona-fide CCs and allows us to rank confi- 
dence in membership according to membership probabilities. 
Our approach is 8 dimensional in the sense that 3 spatial, 3 
kinematic, and two population parameters (iron abundance 
and age) are used as membership constraints. Both data in- 
homogeneity and incompleteness are critical limitations to 
this work, and are acknowledged in the relevant sections. We 
describe our analysis in Sec. [2] 

For the first time, we systematically search for cluster 
members among Cepheid candidates from surveys such as 
AS AS, NSVS, ROTSE, and also from the suspected vari- 
ables in the General Catalog of Variable Stars. Most data 
employed to do so are taken from published catalogs or other 
literature. However, we also perform radial velocity observa- 
tions of Cepheids on both hemispheres and determine sys- 
temic velocities, Vj. To improve sensitivity to binarity, lit- 
erature RVs are added to the new observations. The data 
compilation is described in Sec. [3] 

The bona-fide CCs identified by our analysis are pre- 
sented in Sec. 14. II and certain cluster- Cepheid combinations 
(Combos) are discussed in some more detail. A full table 
containing the cluster Cepheids candidates identified in this 
work is provided in digital form through the CD^. Com- 
bos that deserve observational follow-up are identified in the 
text. Particular attention is given to Combos previously dis- 
cussed in the literature. Discussions of additional Combos 
can be found in the online appendix. In Sec. l4.2l we then em- 
ploy our bona-fide CC sample in a calibration of the Galac- 
tic Cepheid PLR. The method and results are discussed in 
Sec.[5l which is followed by the conclusion in Sec.|6l 



2 MEMBERSHIP ANALYSIS 

Our all-sky census is structured as shown in Fig.[T] First, 
lists of known open clusters and known Cepheid candidates 
are compiled, see Sec.[3] for details. Second, the two lists 
are cross-matched positionally in a many-to-many relation- 
ship so that we investigate a given Cepheid's membership 
in multiple different open clusters, and a given open cluster 
can potentially host multiple Cepheids. The correct classi- 
fication of cross-matched Cepheid candidates is verified by 
considering light curves, and spectra. Misclassified objects 
are removed from the Cepheids sample. Third, membership 
probabilities are calculated based on all available member- 
ship constraints. These last two points are described in the 
present section. 

Membership probabilities are cal culated follo wing 
Bayes' theorem that can be formulated as (|javnej2003l . M): 

P{A\B) = oc P{B\A) X PiA). (1) 

The posterior probability P{A\B) ('membership probabil- 
ity') is proportional to the product of likelihood, P{B\A), 
and prior, P{A). P{B\A) represents the conditional proba- 
bility of observing the data under the hypothesis of member- 
ship, and P{A) quantifies the degree of initial belief in mem- 
bership. The normalization term P{B), of which we possess 
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Figure 1. Schematic view of membership analysis. Rectangu- 
lar boxes represent data sets used, green rounded boxes indicate 
actions. Cepheids are cross-matched with open clusters to form 
Combos. Only open clusters with consistently defined core and 
limiting radii are considered. Data from the literature and new 
observations are combined for each Combo. Cepheid classification 
is verified based on the data compiled (light curves, spectra). Pri- 
ors, P{A), and likelihoods, P(B\A), are calculated separately and 
joined as membership probabilities, P(A\B). 



no knowledge, is the probability to observe the data. We 
define P{A) in Eqs.|3]&|land P{B\A) in Eq.[Z] below. 

2.1 Prior Estimation and Positional Cross- match 

2.1.1 Positional Cross-match 

On-sky proximity is a necessary, but insufficient criterion for 
membership. Intuitively, if no other information is available, 
one might tentatively assume membership for a Cepheid 
that falls within the core radius of a potential host cluster. 

Therefore, our census starts with a positional cross- 
match that aims to identify all combinations of cluster- 
Cepheid pairs that lie sufficiently close on the sky to warrant 
a membership probability calculation (Combos). The cross- 
match itself is straightforward: if the separation between a 
cluster's center coordinates and the Cepheid's coordinates 
is less or equal to 5 limiting cluster radijf], we include the 
Combo in our analysis. Using this proximity criterion, we 
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^ This cut-ofF radius was adopted to include possible members of 
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cross-match 583 different open clusters (from 1658 initially 
compiled) with 1028 Cepheids (from 1844 initially compiled) 
and obtain 2570 Combos that we investigate for member- 
ship. 

The initial cross-match is purely positional, and the ma- 
jority of Combos studied are non-members. Our analysis in- 
tends to weed out this majority and to indicate to us the 
good candidates through a high membership probability. 

2.1.2 The Prior 

We define the prior, P{A), using the on-sky separatior[f| be- 
tween cluster center and Cepheid, weighted by the cluster 
radius, i.e. its apparent size on the sky. 

The radius of an open cluster is typically determined 
by fitting an exponential radial density profile to a stellar 
over-density on the sk y, an approa ch originally developed 
for globular clusters bv iKind (| 19621 ). The method relies on 
the assumption that two separate distributions are seen: a 
constant field distribution and one that is attributed to the 
cluster. 

Various ways to define cluster radii can be found in 
the literature. Among these are the 'core radius' (most 
stars belong to cluster), rc, and the 'limiting radius' 
for the cluster halo (str o ng fiel d star contamination), 
'"lim, see iKharchenko et al.l l|2005bl . from hereon: K05) and 
iBukowiecki et al. I l|201ll . from hereon: Bll). 

Intuitively, the probability o f membership i s relat ed to 
separation and cluster radius, cf. ISanchez et all (|2010l ). Let 
us therefore define the quantity x as: 

(2) 



2nim - Tc 

where r denotes separation, x is negative, if the Cepheid lies 
within the cluster's (projected) core and becomes unity at a 
separation equal to twice the limiting radius. We define our 
prior, P{A), so that (no other constraints considered) mem- 
bership is assumed when the Cepheid lies within the cluster's 
core, i.e. a; < 0. Outside rc, inspired by radial density pro- 
files of star clusters, we let the prior fall off exponentially 
and define it to reach O.lpercent = 10~^ at a:: = 1. Hence, 



P{A){x < 0) 
P{A){x ^ 0) 



1 

10" 



(3) 
(4) 



Figure[2] serves to illustrate this definition. The prior thus 
carries the 2-dimensional information of separation and clus- 
ter radius, and thereby takes into account how concentrated 
a cluster is on the sky assuming circularly distributed mem- 
ber stars. 
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Figure 2. Illustration of the adopted prior, P{A), as a func- 
tion of separation normalized to a cluster's radii, expressed by 
the quantity x, cf. Eq.[2l If a; (Cepheid within core radius, 
rc): ^"(^4) = 1. Outside the core, P(A) decreases exponentially, 
inspired by radial density profiles of stars clusters. We adopt 
P{A) = 0.001 at r = 2riim. 



approach w as inspired by t h e Hip parcos astrometry-based 
studies by iRobichon et al.l (|l999l ) and iBaumgardt et al.l 
(200^. We extend it here to take into account up to 6 dimen- 
sions using parallax, tt, radial velocity (RV), proper motion, 
and fis, iron abundance, [Fe/H], and age (open clusters 
assumed to be co-eval), weighting all constraints equally. 

Assuming that a given Cepheid was not used to deter- 
mine a cluster's (mean) parameters, we can calculate the 
quantity 

c = x'^T,~^x, (5) 

where x denotes the vector containing as elements the dif- 
ferences between the (mean) cluster and Cepheid quantities: 

X = (tTCI - TTCcp, («r,Cl) " «7,Ccp, . . ■) . (6) 

Let Cci be the covariance matrix of the cluster and Ccep 
that of the Cepheid. Let S then denote the sum of the two 
and its inverse. Since the data employed in this cal- 

culation comes from many different sources, no knowledge 
of correlations between the different parameters is avail- 
able. We thus make the assumption of independent mea- 
surements, which results in diagonal covariance matrices 
containing only parameter variances. Possible correlations 
between Cepheid and Cluster parameters are thus assumed 
to be negligible. We consider this justified, since we pos- 
sess no knowledge of the extent of such correlations and 
assume that Cepheids were not used in the determination of 
cluster mean values. This formulation furthermore implicitly 
assumes normally (Gaussian) distributed errors. 



2.2 The Likelihood P(B|A) 

The likelihood, P{B\A), is computed as a hypothesis test. 
It estimates the probability that the observed data is con- 
sistent with the null hypothesis of (true) membership. This 



cluster halos in the analysis, inspired by the well-known case of 
SZ Tau in NGC 1647. 

^ We avoid the term 'distance' when referring to the on-sky sep- 
aration (in arcmin) in order to prevent confusion with radial dis- 
tance (in pc) 



XjVj^f , where A^dof is the number of degrees of freedom equal 
to the length of vector x, ranging from 1 to 6. c thus de- 
pends on the number of membership constraints considered 
(the on-sky position is used in the prior). In cases where no 
membership constraints are available, i.e. A^dof ~ 0, we set 
PiB\A) = 1. 

P{B\A) is obtained by calculating unity minus the p- 
value of c, p(c): 



PiB\A) 



p(c). 



(7) 



Since the distribution (and therefore the p-value com- 



4 R.I. Anderson et al. 



puted) is very sensitive to A^'dof for small A^dof, P{B\A) nat- 
urally contains information on the number of membership 
constraints employed. 

Of course, we cannot prove the null hypothesis, only 
exclude it. However, by including the greatest number of 
the most stringent membership constraints possible, this 
method very effectively filters out non-members. The re- 
maining candidates can therefore be considered bona-fide 
members, provided the constraints taken into account are 
sufficiently strong. 

The filtering effectiveness of the likelihood strongly de- 
pends on the uncertainties adopted for the constraining 
quantities: the larger the error, the weaker the constraint. 
Conversely, the smaller the error, the more important be- 
come systematic differences between quantities measured or 
inferred through different techniques. Obtaining reasonable 
estimates of the external uncertainties is of paramount im- 
portance to the success of this work, since the data con- 
sidered is inhomogeneous and listed uncertainties typically 
provide formal errors or estimates of precision. 

For certain quantities, we therefore adopt increased er- 
ror budgets that we motivate and detail in the following 
sections. Care is taken to avoid too large or too small error 
budgets, and to ensure that likelihood remains an effective 
membership criterion. 



3 DATA USED TO COMPUTE LIKELIHOODS 

In this section, we describe how we compile the data used 
for our analysis. The constraints employed are: on-sky sepa- 
ration, parallax, proper motion, radial velocity, and the pop- 
ulation parameters iron abundance (as a proxy for metallic- 
ity) and age. Most data considered originates from published 
literature and catalogs. Ifowever, we also include radial ve- 
locity (RV) data from an extensive, year-round observation 
program carried out on both hemispheres. Some details on 
this program are provided in Sec. 13. 2. 3. II A full description, 
however, is out of scope for this work and will be published 
separately. 

Very often, data on a given membership constraint can 
be found in different references. In such cases, a choice of 
which reference to prefer over the other ones has to be made. 
In each of the following subsections, the references men- 
tioned first are the ones preferentially adopted. This section 
is divided into two parts: Sec . 13 . II dedicated to open clusters, 
and Sec. l3.2l to Cepheids. 



3.1 Open Cluster Data 

The definition of P{A) in Eq.|4]sets an important constraint 
on the literature on open clusters to be considered, since 
both a core and limiting radius have to be defined and 
both quantities have to be comparable in nature between 
different catalogs. We found three extensive studies that in- 
clude large numbers of op en clusters and sa tisfy this re- 
quirement: K05, Bll, and iKharchenko et all (^012) (from 
hereon: K12). The study bv lFroebrich et all l|20o'7Fr was not 
included here, since the false positive rate may exceed 50 per 
cent (see also ICamargo et al.lboiol ). and since we noticed 
a suspicious correlation between Vc and rum. Note, however, 
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Figure 4. Left panel: comparison between ri in K12 and the 
core radius, Rco, in K05; right panel: same for r2 in K12 and the 
limiting radius, R^-i, in K05. The radii are comparable. 



that some clusters list ed in K12 were originally identified by 
iFroebrich et"al1 (I2007D . 

Since many clusters in K12 were also studied by K05, 
we compare the three radii defined in K12 with the core 
and limiting radii in K05 and notice that the limiting ra- 
dius in K05, Rci, corresponds well to r2 in K12 (though r2 
tends to be smaller), while ri in K12 is rather similar to 
the core radius in K05, Rco- Nevertheless, a fair amount of 
scatter exists between both studies, see Fig. [J] We consider 
K12 an update (and extension) of K05 and therefore prefer 
the newer cluster parameters over the older ones. 

We previously compared radii given i n K05 and Bll 
for th e clusters common to both works in [Anderson et al.l 
(|201^. Rather large scatter is present (more than a factor 
of 2 for an appreciable fraction) and illustrates that clus- 
ter radii are subject to significant uncertainty. However, the 
radii from both studies follow the same trend and we there- 
fore consider them comparable for our purpose, although 
K05 is b ased on optical and Bll on near-infrared (NIR) 
2MASS (jCutri et al.ll2003l ) photometry. 

Given the sometimes rather large difference between 
cluster radii mentioned in the literature, we adopt a 'per- 
missive' scheme that gives preference to the study giving the 
largest limiting radius for the cluster in question in order to 
take into account the most interesting cases (in number and 
qualitatively speaking) and thereby bias ourselves towards 
higher P{A). We therefore strongly rely on the remaining 
membership constraints that define the likelihood of filter- 
ing out chance alignments. 

We then ad d to our compilation data from the 
iDias et al.l (|2002l ) catalo£| (from hereon: D02), the ex- 
tended Hippajcos c atalo g (from now on: XHIP) by 
lAnderson fc Francis! (I2OI2I) and the o pen cluster parame- 
ters catalog bv Iciushkova et al.l (|2010l ). Following D02, we 
remove the clusters Haffner 3, 5, 25 and H ogg 3 from our sam - 
ple, since they are probably not clusters (|Piatti et al.ll201ll ). 
The final list contains 1658 entries, 1296 of which are open 
clusters that are also listed in D02, 7 of which are OB as- 
sociations from K05, and 355 of which were newly added 
by including K12. 583 open clusters are cross-matched with 
Cepheids and evaluated as potential Cepheid hosts. 

Unf ortunate l y, hal f of the 24 Cepheid-host clusters men- 
tioned in [Turned I2OI0I ) are not present in the catalogs used 
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to compile our list of cluster^. We therefore cannot investi- 
gate these particular cases here. 

FigureO shows the distribution of clusters (blue solid 
circles, scale with limiting radius) and Cepheids (green open 
squares) in Galactic coordinates. Clusters closely trace the 
disk, and no obvious gaps are present in our all-sky census. 

As mentioned in Sec. l2.1l we adopt a 'permissive' 
scheme regarding the open cluster radii. We consistently ap- 
ply this logic and always adopt the center coordinates and 
radii of the reference giving the largest limiting radius. For 
consistency, the parameters distance, reddening, age, etc. 
are preferentially taken from the same source. Distances in 
[pc] (and their uncertainties) are converted to parallaxes in 
[mas] through Eqs. [T2]and[ni see Sec. imi 

In the following subsections, we outline how we con- 
struct the open cluster data set for proper motion, radial 
velocity, iron abundance, and age. 

3.1.1 Mean Proper Motion 

Mean cluster proper motions, p-a.ci ^^'^ fis,ci, are taken 
from K05, K12, and D02. 

The uncertainties on mean proper motion listed in these 
references are typically calculated either as intrinsic disper- 
sions (for clusters closer than approx. 400pc, K12), or as 
standard mean errors, i.e. the error decreases as y/N, — 1, 
where A^, is the number of stars considered members, cf. 

The quoted uncertainties on the cluster mean are thus 
much smaller than the uncertainty on an individual cluster 
star's measurement. For example, in K12, the typical mean 
proper motion error is 0.4 masyr"^. 

For the majority of Cepheids, however, the uncertain- 
ties on proper motion are much larger, and many have been 
obtained from different data sets, using different techniques. 
Therefore, to ensure comparability of inhomogeneous data 
and to reduce our sensitivity to offsets in zero-points due 
to data-related specificities such as reduction techniques, 

5 The missing clusters are: ADS 1477, ADS 5742, AnonVulOB, 
CasR2, Dolidze45, Lynga6, NGC 6649, Ruprecht 175, Trum- 
pler35, Turner 2, VelOB5, VulOBl 

S See under 'version 2.3 (25/abr/2005)' in file: 
http : //www. astro . iag.usp .br/'^wilton/whatsnew. txt 



we adopt a more generous error budget for /i* qi and fls,ci 
that resembles the uncertainty of an individual cluster star's 
proper motion. This is done by multiplying the uncertainty 
listed by the factor VA* — 1 and thus slightly reduces the 
weight of proper motion as a membership constraint. Em- 
pirically, we are confident that this is justified, since proper 
motions of Cepheids typically barely exceed their uncertain- 
ties, and care should be taken not to over-interpret their 
accuracy. 



3.1.2 Mean Radial Velocity 

Average radial velocities for open clusters were taken from 
K05, K12, the extended Hipparcos compilation (XHIP), and 
D02. We adopted the XHIP values where available, and oth- 
erwise adopt the mean cluster velocity based on the largest 
number of starfl 

Qualitative differences exist in the uncertainties listed 
for mean cluster RVs. For some well-studied clusters, the un- 
certainty given is an estimate of the intrinsic RV dispersion. 
The majority of cluster RVs, however, were estimated from 
only a few stars (about half on two stars or less, cf. Fig.[5l 
and are therefore subject to systematic uncertainties due to 
implicit membership assumptions, for instance. In addition, 
unseen binary companions and instrumental zero-point off- 
sets can introduce systematic uncertainties at the level of a 
few kms~^. 

We therefore adopt 2kms^^ as a minimum uncertainty 
of the mean cluster velocity. If no uncertainty estimate is 
given, we adopt a{RVci) = lOkms"'^ / Nrv as a typi- 
cal uncertainty on the mean cluster RV, where Nrv is the 
number of stars used to determine the mean cluster RV. 



3.1.3 Iron Abundance 

We adopt iron abundances compiled in D02, including the 
uncertainties given. The mean uncertainty among the clus- 
ters compiled is 0.08 dex. 



With the exception that the values in K12 replaced the values 
published in K05 
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Figure 5. Cumulative fraction of clusters for which a given num- 
ber of stars, N^y, was used to determine the average RV. About 
half of average cluster RVs are based on a couple of stars. 



3.1.4 Cluster Age 

Ages were available for most clusters, since they are often 
determined simultaneously with the distance via isochrone- 
fitting. Such was the case in Bll, for instance, or for many 
clusters in K05. Although a model-dependent parameter, 
age does provide a valid constraint for membership, reflect- 
ing evolutionary considerations that are empirically vali- 
dated. Quantifying an uncertainty for age as a parameter, 
however, is rather difficult. 

Younger clusters exhibit a Main Sequence turn-off at 
higher stellar masses than older clusters. As a consequence of 
the Initial Mass Function, a younger cluster's turn-off point 
tends to be less populated than that of an older cluster. It 
therefore follows that age estimates tend to become more 
accurate with age, since the cluster's turn-off point tends to 
be defined more clearly against the field and therefore better 
constrains an isochrone fit. 

Figure|6] corroborates the above reasoning by showing 
cluster ages against their uncertainties as given in K12. We 
thus estimate an upper limit on the uncertainty of cluster 
age as the red dashed line in Fig.|6l which is: 



cr(logaci) ^ 0.3 - 0.067 (logaci - 7.0) 



3.2 Cepheid Data 



(8) 



Cepheid candidates were compiled from the January 
2012 version of t he General Catalog of Variable Stars 
ijSamus et al.ll20ll from hereon: GCVS) and the May 2012 
version of the AAVSO Variable Star Index (from hereon: 
VSXfl From GCVS and VSX, we import the variability 
types CEP, CE P(B), DCEP DCEPS; fr om VSX, we in- 
clude the AS AS iPoimanskil (|l997l . \200i ): iPoimanski et all 
l|2005ll Cepheid candidates classified as DCEP-FU or DCEP- 
FO. This li st also contains C epheid ca ndidates found 
by R OTSE l|Akerlof et al.ll2000l ) or NSVS (|Wozniak et all 
I2OO4'), as well as the ones in the suspected variables catalog 
(,Kukarkin & Kholopov 1982). 

This starting point contains an unknown, but probably 
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Figure 6. Age uncertainty, (T(logaci), as function of cluster age, 
logaci, given by K12. Older clusters have more precisely esti- 
mated ages. We adopt as error budget for clusters without stated 
age uncertainties an upper limit to this proportionality indicated 
by the red dashed line, cf. Eq.[8] 



high, fraction of non-Cepheids. Contamination due to type- 
II Cepheids, or Cepheids belonging to the Magellanic Clouds 
are removed from the sample during the cross-match with 
the clusters (trace the disk). To further reduce contamina- 
tion, we visually inspect all ASAS-3 V-band fight curves of 
Cepheid candidates with ASAS identifiers. 

Radial pulsation and color variations during the pulsa- 
tion are defining characteristics of Cepheids. We thus use the 
spectra obtained for radial velocity observations described 
in Sec. l3.2.3.1l to verify classification. A total of 146 ASAS 
Cepheid candidates and 18 others are thus rejected from the 
Cepheid sample, resulting in a final list of 1844 Cepheid can- 
didates, 1028 of which are cross- matched with open clusters. 

The cleaned sample of Cepheids cross-matched with 
clusters was appended with literature data from many 
sources, and references are given in the text. Among the 
most relevant references are: 



• The lFernie et al.l (Il995l) DDO Cep heid databasfH 

• The Klagvivik &: SzabadosI (|2009l 'l Cepheid database 
(KS09) 

• T he ASAS Catalog of Variable Stars (jPoimanski et akl 
l2005l . ACVS) and associated photo metry 

• The new Hipparcos reduction l|van LeeuwenllioOTi ) 

• The extended Hipparcos compilation 

d Anderson fc Francis! [201^ . XHIP) 

• The ASCC-2 5 cata log (jKharchenkdl200ll ) updated by 
iKharchenko et al.l (120071 ') 

• The PPMXL catalog |Roeser et al."2010|) 

• The 2MASS catalog (Cutri ct al. 2003) 

• The Cepheid photom etry obtained by Berdnikov et al.l 
(|2000l ): lBerdnikovl (|2008^ 

• The McMaster Cepheid photometry and radial velocity 
data archive maintained by Doug WelclF°l 

• The radial velocity data, see Sec. l3.2.3l 



http : / /www . aavso . org/vsx/ 



^ http : //www. astro .utoronto . ca/DDO/research/Cepheids/ 
http : //crocus .physics .mcmaster . ca/Cepheid/ 
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3.2.1 Cepheid Parallaxes 

Parallax, tt, is a key membership constraint, since cluster 
membership is virtually guaranteed if a Cepheid occupies 
the same space volume as a cluster. We combine parallax es- 
timations from different sources, favoring PLR-independent 
determinations. 

Parallax in [mas] is given preference over distance in 
[pc] here, since the uncertainty, u-^, is normally distributed, 
in contrast to the error in distance. This is important, since 
the computation of likelihoods by Eq. [7] assumes Gaussian 
uncertainties. 

We compile parallaxes in the following order of prefer- 
ence: i) Hipparcos absolut e parallaxes and asso ciated errors 
from the new reduction by Ivan Leeuwenl l|2007l . 5 Cepheids), 
so lo ng as a-n/T^ ^ . 1 and tt > ; ii) HST- based parallaxes 
from [Benedict et al.l (|2007, 8 Cepheids ) : iii) parallaxes com- 



, . , p a 

puted from distances in Istorm et all (|201ll . 67 Cepheids) 
including their uncertainties; iv) parallaxes calculated from 
PLR-based distances (624 Cepheids), see below. 

PLR-based parallaxes of fundamental-mode Cepheids 

are calculat ed from distances computed following 

[Turner et all (|201(]| ). Our choice of P-L relation was 
motivated mainly by the considerations that i) V-band 
magnitudes can be obtained for the largest number of 
Cepheids; ii) the above formulation is calibrated for the 
Galaxy using the most rece nt observational results, includ- 
ing the HST parallax es by 'Benedict et al.l (|2007l ') and the 
cluster Cepheids from [Turn er (2010). 

We thus calculate PLR distances as follows: 



5 log d = {mv) - (Mv) - Ay + 5 , 



(9) 



where (mv) is the apparent mean V-band magnitude, and 
the average absolute V-band magnitude, (Mv), is obtained 
from the pulsation period P via: 



(Mv) = - (1.304 ± 0.065) - (2.786 ± 0.075) log P . 



(10) 



Eq.[TO]is valid only for fundamental- mode pulsators, no dis- 
tances were estimated for overtone pulsators. The total ab- 
sorption, Av, is defined as: 



Av ^Rv E(B - V) , 



(11) 



with Rv = 3.1 the canonical ratio of total to selective ex- 
tinction (reddening law) and E{B — V) the color excess of 
the object, cf. Sec. l3.2.1.2l The parallax is simply: 



1000 

d 



(12) 



with d in [pc]. The parallax uncertainty, cr^, is obtained con- 
sidering the error budget on the distance, ad'. 

1000 



d2 



■ aa [mas] . 



(13) 



Thus, to estimate a Cepheid's parallax, knowledge of 
the PLR, P, (mv), and Ay is required. Periods are usually 
available in the GCVS or the VSX, whereas average magni- 
tudes and color excesses of many of the newer Cepheid can- 
didates are not available in the literature. However, E(B-V) 
can be estimated from combined NIR and optical data. The 
following paragraphs describe in detail how these quantities 
are compiled. 



120 
100 
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(^KS09,f) 



0.00 



(Akso9.b) = 0-03 



Figure 7. Histogram of differences in mean magnitudes rela- 
tive to the KS09 values, Akso9- Blue slim bars show Akso9,f 
computed using Fernie values (123 Ceph eids); Red broad bars 
show Akso9 B computed using data from [Berdnikov et al.lBoOOl 
(127 Cepheids). Mean differences, (AkS09)i ^rid dispersions, 
o-(Akso9). given in [mag]. 




Figure 8. Histogram of differences in mean magnitudes rela- 
tive to Hipparcos median V-band, Ajjip. Red slim bars show 
Ahip,kS09 computed using KS09 ( 104 Cepheids). Blue b road 
bars show Ajjip b computed using [Berdnikov et al.l I2OO0I (198 
Cepheids). Mean differences, (Ajjip), and dispersions, (t(Ahip), 
given in [mag]. 




Figure 9. Histogram of differences in mean V-band magnitude 
computed from ASAS light curves for 154 Cepheids relative to 
reference values from Fernie and KS09, AasaSi cf. also Fig.[7l 
Mean difference, {Aasas)i s^nd dispersion, (t(Aasas)i given in 
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3.2.1.1 Mean Magnitude, (mv) We compile mean 
V-band magnitudes, (mv), from multiple references. Dif- 
ferent methods of determining mean magnitudes exist, 
and the photometry employed is inhomogeneous, forcing 
us to adopt a zero-point for mean magnitudes compiled. 
In Fig.[7| we therefor e com pare mean magnitudes from 



Klagyivik fc Szabadoj (|2009l . from hereon: KS09) with the 



Fernie et al.l ((1995) datab ase's magnitude-based m eans and 



the intensity- means from iBerdnikov et al.l l|2000h . For the 
Cepheids common to both studies, KS09 and the Fernie 
magnitudes show excellent agreement. We therefore adopt 
the following order of preference for compiling mean magni- 
tudes. 

First, we adopt {mv) values from KS09 with a fixed er- 
ror budget of 0.03 mag, since the study carefully investigates 
amplitudes with a special focus on binarity. 

Second, we adopt the magnitude-based means from the 
Fernie database with uncertainties calculated as the differ- 
ence between intensity- and magnitude-mean magnitudes, 
with a minimum error of 0.03 mag. 

Third, we include IBerdnikov et all (|2000l ) mean mag- 
nitudes.. As seen in Fig.[71 these {mv) values are system- 
atically smaller (brighter) by approx. 0.03 mag than KS09. 
This discrepancy is most likely due to different ways of 
determining the mean. We remove this offset from the 
IBerdnikov et ak 

1 (1200(f) values for internal consistency and 
adopt 0.03 mag as error budget for these values, identical to 
KS09. 

Fourth, we empl oy median V-b and magnitudes from 
the Hipparcos catalog (|Perrvman fc E SA 1997, obtained via 
XlflP) for 8 Cepheids. The median V-band magnitudes 
derived from Hipparcos magnitudes can differ significantly 
from mean magnitudes listed in other references, cf. Fig. [8] 
Usually, this is due to contamination due to a nearby star 
within the instantaneous field of view. As error budget for 
the Hipparcos median V-magnitudes, we adopt the stan- 
dard deviation computed from the red histogram in Fig.|8l 
i.e. 0.110 mag. We note that we could find no dependence 
on period or number of transits for this dispersion. 

Fifth, we adopt average apparent V-band magnitudes 
that we determine from ASAS-3 light curves. To this end, 
we fit Fourier Series (same procedure as described for radial 
velocities in Sec. l3.2.3.3|) to the phased light curves and use 
the constant term as the average, (mv). Figure [5] shows a 
histogram of Aasas, the differences between the computed 
ASAS-based (mv) and Fernie or KS09. We remove the offset 
of —0.01 mag from the ASAS mean magnitudes and adopt 
the dispersion of 0.10 mag computed as the error budget. 

The large dispersion, a (Aasas), in Fig.[9]probably orig- 
inates from contamination due to nearby stars. To illustrate 
this. Figure [TT] shows phase-folded ASAS-3 V-band light 
curves of two Cepheids, CYCar (left) and BMPup (right). 
Our mean magnitude agrees well with the literature value 
for CYCar. For BMPup however, a systematic difference 
of approximately 0.144 mag is evident, although the light 
curve appears to be clean otherwise. Inspection of a DSS 
images, however, reveals that contamination from a nearby 
companion is likely. Out of 154 Cepheids for which the ASAS 
light V-band curves were inspected, 20 differed by more than 
0.1 mag from the reference value, and 28 agreed to within 
0.01 mag. 

If no mean magnitude is obtained from any of the 



g 60 




Figure 10. Fractional amplitudes, fa, computed from Fernie 
mean, and GCVS maximum and minimum magnitudes, see 
Eg .[Til 



above sources, we perform a (rough) estimate of (mv) based 
on the information provided in the GCVS and the VSX, 
using the magnitude at maximum brightness, miuv, and 
the amplitude, amp^, of the V-band light curve, amp^ is 
either provided directly by the catalogs, or calculated as 
the difference between minimum and maximum brightness, 
amp^ = maxv — minv- 

Since Cepheid light curves are skewed, their mean mag- 
nitudes do not necessarily lie at half the amplitude. We 
therefore estimate the typical fractional amplitude at mean 
brightness, {fa), to compute (mv) = minv + (/a)ampy. 
Figure [10] shows a histogram of fa computed using mean 
magnitudes listed in the Fernie database, (mv.p), and am- 
plitudes, amp^, from the catalogs. We find 



{fa) = median 



(my^p) — maxv 
minv — maxv 



0.54 ±0.079. 



(14) 



We derive an uncertainty on (mv) thus obtained using 
the uncertainty on {fa), an estimated error on the ampli- 
tude, and a prescribed error on the magnitude at maximum 
brightness (we adopt 0.1 mag for 12th magnitude and be- 
low, and increase linearly to 0.5 mag at 20th magnitude). 
The resulting mean error on {mv) is 0.27 mag. 

This estimation works reasonably well, although there 
exist obvious limitations, such as inhomogeneity of pass- 
bands, accuracy of the upper and lower limits, the applica- 
bility of the above ratio for a given Cepheid. Nevertheless, it 
does provide access to rough estimates of (mv) for Cepheids 
with little available information. 

Given the many different ways in which average magni- 
tudes were estimated, we keep track of the type of estimation 
to ensure traceability of any potential issues. 

3.2.1.2 Color excess, E(B — V) The principal refer- 
ence adopted for color excess is iKovtvukh et all (120081 ) from 
which we also adopt stated uncertainties, cr(E(B — V)). 
Where no o (E(B — V)) were listed, we adopt an error bud- 
get of 0.05 mag, which is the mean of the giv en uncertainties. 

F or Cepheids without measurements in lKovtvukh et al] 
(|2008l ). we adopt E(B-V) from the Fernie database and an 
error budget of 0.05 mag, unless the standard error from 
multiple reddening estimations was given. 
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Figure 11. Phase- folded ASAS-3 V-band light curves of CYCar (left) and BM Pup (right). Julian date of observation is color-coded, 
increasing from red to blue. Horizontal lines indicate reference average magnitude (red dashed, KS09) and constant term of the fitted 
Fourier Series (blue solid). For CYCar, the two are in excellent agreement. BMPup has a bright neighbor that contaminates the aperture 
used to measure its flux, leading to an underestimated (too bright) mean magnitude. 
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Figure 12. Color excess from the literature against 2MASS- 
based estimate. Red solid line indicates the diagonal, cyan colored 
dashed line represents weighed least squares fit. RMS around ei- 
ther line: 0.10 mag. 



F urther E(B-V) values were adopted from lFougue et al] 
l|200'iD together with the stated uncertainties. 

Color excesses for other Cepheids are estimated using 
single-epoch 2MAS S JCutri et aLlbOOSl) J -band magnitudes, 
mj (JD), following iMaiaess et al.l (|2008h . The method re- 
quires knowledge of the pulsation period, P, and the average 
J-band magnitude, (mj). To obtain the latter, we employ 
their Eq. 5 that reads: 



(mj) ~ m.j(JD)- 



\mv(4,,) - maxvl 



amp^ 



-0.5 



•O.STampv- , (15) 



where (j>J denotes the pulsation phase of the J-band mea- 



surement, and m 



V(0.7) 



is the V-band magnitude at that 



phase. E(B — V) can then be estimated by 



E{B-V) = -0.270 log P-f 0.415 ((mv) - (mj))-0.255 .(16) 

Wherever possible, mv (,),,) was obtained from the ASAS 
light curve. If this is impossible, we assume a sinusoidal light 
curve with the given mean magnitude and (semi-)amphtude. 

Uncertainties or changes in pulsation period can sig- 
nificantly impact the phase calculated for the single-epoch 
2MASS measurement, (j>j. We therefore optimize Cepheid 
ephemerides for which ASAS data were available. To do so, 
we compute a grid (at fixed periods) of Fourier Series fits 
around the period provided in the ACVS and retain the so- 
lution with the minimum root mean square. Epochs are op- 
timized by simply shifting the phase-folded curve. Figure [T3l 
illustrates this step for the overtone Cepheid QZNor. We 
then take care to employ the most recently determined pul- 
sation ephemerides available and estimate reddening uncer- 
tainties using error propagation for the quantities involved. 

We note that this approach may be subject to multi- 
ple issues such as: i) the unknown shape of the light curve; 
ii) the applicability of Eq.[T6l iii) period changes that im- 
pact mv(^j); iii) the approximate form of the relationship 
in Eg . 1151 We therefore compare the 2MASS-based color ex- 
cesses to the reference values, see Fig. ll2l where the result 
of a weighted least squares fit is indicated by a dashed cyan 
line and does not differ much from the diagonal indicated 
by a solid red line. Despite considerable dispersion, the cor- 
respondence is clear and the results are promising (RMS of 
0.1 mag). 



3.2.2 Proper Motions 

Cepheid proper motions are taken from the following sources 
in order of preference: 

(i) Hipparcos pro per motions from the new reduction by 
Ivan LeeuwenI (|2007l ) 

(ii) The PPMXL catalog bv lRoeser et~aLl l|2010l) 



3.2.3 Systemic Radial Velocities 

3.2.3.1 New Observations In order to extend the 
number of Cepheids with known systemic radial velocities. 
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Figure 13. ASAS-3 V-band phase-folded light curve of QZNor. Color-code from red to blue indicates increasing Julian date of obser- 
vation. Left panel: pulsation period and epoch of maximum light from the ACVS. Right panel: optimized period and epoch used. 



f-y, wj"! carried out observations between November 2010 
and July 2012 using the fiber-fed high-resol ution echelle 
spectrographs CORALIE (|Queloz et al.l'200ll see also the 
instrumental upgrades described in fScgransa n et al.l |2010| . 
R ~ 60000) at t he 1.2 m Euler teles cope at La Silla, Chile, 
and HERMES (|Raskin et al.l l201ll . R ~ 80000) at the 
identically-built Mercator telescope on La Palma. In total, 
we observed 103 Cepheids with CORALIE and 63 with 
HERMES. 18 Cepheids were observed with both instru- 
ments, i.e. from both hemispheres. For 85 of these Cepheids, 
no radial velocity (RV) data are available in the literature. 

Efficient reduction pipelines exist for both instruments 
that include pre- and overscan bias correction, cosmic re- 
moval, as well as flatfielding using Halogen lamps and back- 
ground modelization. ThAr lamps are used for the wave- 
length calibration. 

The RVs are c omputed via the cros s-correlation tech- 
nique described in iBaranne et al.l (|l996l ). We use numeri- 
cal masks designed for solar-like stars (optimized for spec- 
tral type G2) for all cross-correlations. Both instruments are 
very stable and yield very high precisio n RVs of ~ 10ms~^ 
jQueloz et al.ll2000l : [Raskin et aLlbOllI ). The measurement 
uncertainty is therefore not limited by the instrumental pre- 
cision, but by line asymmetries due to pulsation. A detailed 
investigation of these effects is out of scope for this paper 
and will be presented in a future publication. The typi- 
cal uncertainty on individual measurements is thus at the 
100 — 300 ms~^ level, depending on the star and pulsation 
phase. 

3.2.3.2 Literature Data In addition to the previously 
unpublished radial velocities described in Sec. l3.2.3.1l we 
employ literature data from many references to determine 
systemic velocities, v^, see Sec. 13. 2. 3. 31 The addition of lit- 
erature RVs extends the baseline of our otherwise relatively 
short (1.5 years) observing program, thereby enhancing our 
sensitivity to binarity. For binary Cepheids with published 
orbital solutions, we adopt the literature v-y. 

Aside from the systemic RVs in the Fernie database 
and KS09, we compile RV time series from the follow- 



ing so urces: iLlovd Evanj lll980l): [Cierenl (Il98ll. 1198511 



Coulson et al 



-. Il985h : ICoulson fc Caldwelll 
Barnes et all ll 19871. Il988|); [Cieren et al.l (119891); 
19891): iMetzger et all (1991. 1992, 199|): 



Wilson et al ■ .. . 

Gorvnva et al.l (Il992l'): lEvans fc Welch 



1 1994 1199?! ) 
1 19991): [storm et al 



,1993); Pont ct al 
Bersier et all l[l994l ) : iKisa (.1998. ): .Imbert 
Barnes et al.l (|2005l) : 



2004) 



iPetterson et al.l l|2005l ) : baranowski et all (|2009l ). The 



data for most of these sources dated earlier than 1986 are 
extracted from the McMaster Cepheids databas43- Newer 
data are obtained through VizieBo. 

3.2.3.3 Systemic Radial Velocities, v.^ The systemic 
radial velocity, v^, is obtained by fitting a Fourier series to 
the RVs. We use pulsation period, P, as a fixed parameter, 
since it is known for all Cepheids we observed. 

The basic analytical form applied was a Fourier Series 
with n harmonics and phase (j> is: 



FS„ = + 



E 



fln sin {2nn(j>) + b„ cos (2nn(j>) (17) 



Since the number of data points available varies for each 
star, we do not fix the number of harmonics in this fit. In- 
stead, we iteratively increase the degree of the Fourier se- 
ries until an F-test indicates an overly complex representa- 
tion, i.e. when spurious fit improvement is more likely than 
0.27per cent. For some stars, we therefore use only a simple 
sine function, whereas stars with many measurements are 
fitted using up to five harmonics. We show two examples 
of newly observed RV curves in Figs. [14] and 1151 A full de- 
scription and publication of the new radial velocity data will 
follow in the near future. 

We adopt a fixed error budget of 3kms~^ on Vj. Al- 
though this may overestimate the uncertainties for some 
very good cases, it is intended to account for a range of 
systematic errors, such as unseen binarity, instrumental zero 
point differences and insufficient phase coverage. We are con- 
fident that this error budget is sufficiently large to prevent 
the exclusion of good member candidates, while being suffi- 



The authors are grateful to the observers who contributed in 
this effort; names are given in the acknowledgments. 



maintained by Doug Welch, available 
http: / /crocus .physics .mcmaster . ca/Cepheid/ 
http : / /vizier . u-strasbg . f r/ viz-bin/VizieR 
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Figure 14. RV data recently published in ISzabados et al.ll2012l for first-overtone pulsator GH Car, obtained in southern hemisphere 
with CORALIE as part of our program. Color-code (red to blue) and size-code (larger to smaller) indicate increasing Julian Date of 
observation. Fit of triple-harmonic Fourier Series indicated by black dashed line, indicated by green dashed horizontal line. Top panel: 
phase-folded, Bottom panel: non-folded, showing zooms of the four observing runs during which data were obtained. 
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Figure 15. New RV data for fundamental-mode Cepheid V2340Cyg, obtained in northern hemisphere with HERMES. Color-code (red 
to blue) and size-code (larger to smaller) indicate increasing Julian Date of observation. Fit of triple-harmonic Fourier Series indicated 
by black dashed line. Vj indicated by green dashed horizontal line. Top panel: phase-folded, Bottom panel: non-folded, showing zooms 
of the three observing runs during which data were obtained. 
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ciently small to provide a stringent constraint. If insufficient 
data points render the Fourier fit unsatisfactory, we deter- 
mine a rough estimate of v-, and its error budget by eye. 

3.2.4 Iron Abundance 

We re ly mostly on iron abundances by iLuck fc Lambert 
l|201ll) and complement these with the compilation in KS09 
that made an effort to homogenize measurements from 
many references, notably Andrievsky and collaborators (e.g. 
iLuck fc Andrievskvll2004l ). Unfortunately, the adopted stan- 
dard value of solar iron abundance can vary among refer- 
ences, possibly introducing systematic offsets between dif- 
ferent authors' studies. Furthermore, an estimation of the 
iron abundance in a Cepheid is more complex than in a non- 
pulsating star, since the stellar parameters (e.g. temperature 
and turbulence) vary during the pulsation cycle, and since 
the atmosphere is not static. We therefor e adopt generous 
error budgets, 0.1 dex for the values from lLuck fc Lambert! 
l|201ll) . and 0.15 dex in [Fe/H] for the compiled list in KS09. 

3.2.5 Age 

Cepheid ages can be calculated for first overtone and fun- 
damental mode pulsators using the period-age (PA) rela- 
tions given in .Bono et al., (,2005 ) . For fundamental mode 
Cepheids, we use logt = (8.31 ± 0.08) - (0.67 ± 0.01) log P. 
For overtone pulsators, the relation used is logt — 
(8.08 ± 0.04) - (0.39 ± 0.04) log P. Age error budgets are cal- 
culated from the uncertainties stated for slope and intercept. 



4 RESULTS 

This section presents the results from our census. We start 
with the bona-fide CCs (Sec. l4.l| l. followed by strong can- 
didates that lie far from cluster center and therefore have 
a low prior (Sec. 14. l!^ . For brevity of the main body, we 
defer presentation of inconclusive and rejected Combos that 
have been previously discussed in the literature, as well as 
cases of fore- and background contamination, to appendix 
A. A full table that includes all data used is available online 
through the CDS. 

Each subsection starts with Combos previously dis- 
cussed in the literature. We take care to include relevant 
references. Ifowever, in a field with this much history, it is 
almost inevitable that some works were overlooked. 

As a sanity check of our results, we compare them to 
known CC candidates from the literature. The main refer- 
ences used for this benc hmarking are the reviews bv iFeasj 
(|199 ^ and'Turneil (|2010l . from hereon: TIO). However, there 
are further Cepheids whose cluster membership was consid- 
ered in the literature that will be mentioned where appro- 
priate. As mentioned in Sec l3.1l 12 of the host clusters listed 
in TIO are not present in our cluster sample. For 8 of these, 
the Cepheids are cross-matched and investigated for mem- 
bership in different clusters. However, none of these host 
cluster 'alternatives' was found to be a good match. Among 
these, SWVel has pre viously been in vestigated for member- 
ship in NGC2660 bv iTurner et all (1993) who found the 
same null result. For the Cepheids XCyg, ^Gem, V367Sct, 



Table 1. Membership probabilities of bona-fid e cluster C ep heids. 
Impo rta nt references: aillrwin 1955, b: Kholo povl Il956l . c:lFeast 
1957) . djl^ndagc' 1958, c: Eggcn 1980, f: Turner 1982, g: Walker 
19851 . h: lTurncr 1986 . i: Turner et al. 1992.. j: Ma tthews ct al. 1995 
and references th ere in. k:lTurner et alJIiggTT h lHovle et al]i2CI0a 
m: iAn et aljl200"7l . n: lTuriie'in2010r o: first discussed in this work. 
Column Data lists constraints employed 1: separation (relative 
position, cluster radius), 2:it, 3: RV, 4: ^* , 5: fig, 6: [Fe/H], 7: age. 
P(A\B) = P(A) ■ P{B\A) , see Sec.[2] 



Cepheid 


Cluster 


Ref 


Data 


P{A\B) 


PiB\A) 


SU Cyg 


Turner 9 


k,n 


1-5,7 


0.946 


0.946 


U Sgr 


IC 4725 


a,b,c,m 


all 


0.983 


0.983 


DL Cas 


NGC 129 


b,i 


1-5,7 


0.928 


0.928 


CF Cas 


NGC 7790 


d,j 


1-5,7 


0.752 


0.979 


S Nor 


NGC 6087 


a,b,c,h 


all 


0.710 


0.710 


V340 Nor 


NGC 6067 


g,m,n,l 


all 


0.605 


0.605 


CE CasA 


NGC 7790 


d,j 


1-5,7 


0.540 


0.977 


CE CasB 


NGC 7790 


d,j 


1-5,7 


0.522 


0.963 


QZ Nor 


NGC 6067 


e,g 


all 


0.027 


0.952 


V Cen 


NGC 5662 


f,m,n 


1-5,7 


0.005 


0.923 


V379 Cas 


NGC 129 





1,3,4,5,7 


0.000 


0.923 


SX Car 


ASCC 61 





1,2,4,5,7 


0.001 


0.922 


S Mus 


ASCC 69 





1-5,7 


0.004 


0.879 



and a UMi (Polaris) , no candidate host clusters were iden- 
tified within our cluster sample at separations inferior to 5 
limiting radii. 

In Sec. l4.2l we employ our bona-fide CC sample in a 
calibration of the Cepheid period-luminosity relation. 

4.1 Bona-fide Cluster Cepheids 

Table[T] lists the bona-fide CCs recovered by our analysis in 
descending order of membership probability, P{A\B). Col- 
umn Ref contains the first and most relevant references dis- 
cuss membership for a particular Combo, and column Data 
provides the constraints employed. The last two columns 
contain membership probability and likelihood. 

The horizontal line in Tab.[T]separates the clear from the 
not-quite-so-clear CCs, the difference being that the latter 
are situated at larger separations from the cluster, i.e. have 
low priors. Among these, our work suggests the following 
new bona-fide CCs: SMus in ASCC 69, SXCar in ASCC 61, 
and V379 Cas in NGC 129; see Sec.|4T2]for details. 

4.1.1 Bona-fide CCs with high P{A\B) 

All of the CCs in the top part of Tab. [T] are well-established 
in the literature and are well constrained by the available 
information. They have very high likelihoods and high pri- 
ors, i.e. they are rather obvious members. We therefore refer 
to the original references listed for more information, as well 
as to the data table supplied in electronic forrrF^. 

The case of S Nor in NGC 6087, however, deserves some 
closer inspection, since significantly different cluster RVs 
were found in the literature. 

We note that a large fraction of the Cepheids in Tab.[T] 

Link / catalog number / ? 
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are also binaries, see the binary Cepheid database by L. 
Szabados, and references therein. 

4.1.1.1 SNor and NGC 6087 SNor's members hip in 
NGC 6087 was among the first to ever be s usgested (|lrwirj 
Il955^ and confirmed using ra dial velo c ities (^Feast^ll957^ ■ as 
well as the detailed study bv IXurneil l| 19861 ') based on red- 
dening and distance. 

The values for the cluster's mean RV differ great ly be- 
tween D02 (6kms"^), K05 (-9kms~^), and (Feas^ Il95i 
2.0kms~^), which is important consid ering the Cepheid's 
Vj = 2.53 kms~^ (|Groenewegenl [20081 ). We note that the 
value adopted by D02 is measured on the Cepheid itself 
ijMermilliod et all |2008| ). although without taking orbital 
motion into account, and is tfierc forc not suitable as a 
membership constraint. Since iFeastI ij 19571 ) investigated the 
largest number of stars and specifically targeted this cluster, 
we trust that their 2.0 km s""*^ is the best available estimate 
for the mean velocity of the cluster. 

Small differences exist between the cluster and Cepheid 
values for reddening, metallicity, and age. However, these 
differences barely exceed the combined uncertainties and 
may not be significant. 



4-1.2 Bona-fide CCs with Low Priors 

The following Combos require a short discussion, since we 
obtain high likelihoods and low priors, resulting in low 
mem bership probabilitie s. We note that IRSB parallaxes 
from IStorm et al.l (|201ll ) were employed for all but one 
of the Cepheids discussed here. For SXCar, we adopt a 
fundamental-mode PLR parallax using the E(B-V) value 
from the Fernie database. 

4.1.2.1 QZ Nor and NGC 6067 The overtone pulsator 
QZ Nor is a second Cepheid member to the rich open cluster 
NGC 6067 that also contains V340Nor. Its cluster member- 
ship was first considered by lEggenl (|l983l '). Our likelihood 
of 95 per cent is based on all membership constraints. The 
radius determination for the cluster is rather difficult, since 
it is located in the Norma cloud, cf. the atlas page in K05. 
The low prior may therefore be misleading. Alarmingly, the 
IRSB paraUaxes of QZNor (0.74 ± 0.07 mas) and V340Nor 
(0.58±0.09 mas) are nearly consistent with each other. How- 
ever, their difference barely exceeds the combined uncer- 
tainty, and NGC6067's parallax (0.71 mas) lies between the 
two, consistent with both Cepheids within the error bud- 
gets adopted. It is possible that differential reddening could 
be of importance in resolving this conundrum (higher for 
V340Nor which is closer to cluster center). Furthermore, 
iron abundance and radial velocity of both Cepheids are in 
very good agreement. Proper motion does not exceed the 
uncertainties and is therefore of limited use as a constraint. 
We conclude that the evidence combined in the likelihood is 
in favor of membership, which leads us to consider QZ Nor 
a bona-fide member of NGC 6067. 



http : / /www . konkoly . hu/ CEP/nagytab3 . html 



4.1.2.2 V Cen and NGC 5662 V Gen's available con- 
straints (all but [Fe/H]) are in impressively close agreement 
with those of NGC 5662, yielding a very high likelihood (92 
per cent). As for NGC 6067, the determination of the lim- 
iting radius is complicated by a very strongly populated 
area. Furthermore, the common proper motion clearly shows 
the cluster to be co-moving with the Cepheid. The slightly 
higher proper motion of the Cepheid does not appear to de- 
viate significantly from that of the cluster. We therefore see 
no reason to doubt V Cen's membership in NGC 5662. 

4.1.2.3 SMus and ASCC 69 The membership con- 
straint s considered (all but [Fe/H]) for the binary Cepheid 
S Mus (|Petterson et al.ll2005l ') and ASCC 69 yield a very high 
likelihood of 86 per cent. ASCC 69 is a sparsely populated 
cluster for which K05 list merely twelve 1-a members. The 
radius of the cluster is therefore rather imprecise. Further- 
more, radial velocities are only of limited value as member- 
ship constraints, since the average cluster RV is based on 
only 2 stars. However, proper motion clear indicates that 
cluster and Cepheid are co-moving. We tentatively accept 
S Mus as a cluster member and stress the need for a detailed 
study of its candidate host cluster ASCC 69. 

4.1.2.4 SXCar and ASCC 61 The 4.86 d Cepheid 
SXCar is seen to be co-moving with ASCC 61 in proper 
motion. Unfortunately, no average cluster RV is known. Par- 
allax and age, however, agree very well between cluster and 
Cepheid, lending support to the hypothesis of membership 
with P{B\A) = 95 per cent. An in-depth analysis of the 
cluster, including its mean RV, is of the essence, since the 
cluster is located in a crowded field (K05). We tentatively 
consider SXCar as a member of ASCC 61. 

4.1.2.5 V379 Cas and NGC 129 This high likelihood 
pair at large separation (44' or 3.4 pc at the estimated dis- 
tance to NGC 129) has nearly vanishing proper motion, 
while the both RVs are in excellent agreement and the ages 
are consistent. We obtain a likelihood of 92 per cent. No par- 
allax is computed, since V379 Cas is an overtone pulsator. 
Since NGC 129 has another known member, DL Cas, we can 
compare the pulsational ages of the two Cepheids and find 
both to be consistent within the uncertainties (7.70±0.08 for 
DL Cas and 7.83±0.07 for V379 Cas). Furthermore, the iron 
abundances of both Cepheids are close, as are their RVs (to 
within less than lkms~^). We further note that V379Cas 
and DL Cas have similar reddenin g values with that of 
V379Cas slightly (0.1 mag) higher jKovtvukh et all |2008| . 
both). We therefore tentatively consider V379 Cas a member 
of NGC 129's Halo, pending a better distance estimate and 
additional membership constraints. NGC 129 is thus partic- 
ularly interesting, containing both a fundamental-mode and 
an overtone pulsator, just as NGC 6067. 

4.1.3 Chance Alignment 

Among our bona-fide CCs, 5 are found to lie within Vc, and 4 
a,t rc < R < ri. Within the original cross-match, 21 Combos 
are matched within Vc and 194 within Vc < R < rt. Thus we 
can estimate the rate of chance alignment within the core 
radius to be about four to one (21/5). At separations inferior 
to ri, this ratio increases to approximately 24 to 1 (215/9). 
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Table 2. Membership probabilities calculated for inconclusive 
Combos. Constraints by numbers: l=separation (i.e. position 
&; cluster radius), 2=tt, 3=RV, 4=fi^, 5=fis, 6=[Fe/H], 7=agc. 
P{A\B) = P{A) ■ P{B\A) , see Sec.[2] AH 1 is the abbreviation 
for cluster Aveni- Hunter 1. 



Cepheid 


Cluster 


Constraints 


P{A\B) 


P{B\A) 


BD+47 4225 


AH 1 


1,4,5,7 


0.408 


0.947 


VZ CMa 


Trumpler 18 


1,4-7 


0.000 


0.770 


Y Sgr 


IC 4725 (M 25) 


all 


0.000 


0.702 


GH Car 


Trumpler 18 


1,3,4,5,7 


0.104 


0.574 


BB Sgr 


CoUinder 394 


1-5,7 


0.066 


0.468 


Y Car 


ASCC 60 


1-5,7 


0.337 


0.337 


CG Gas 


Berkeley 58 


1,2,4,5,7 


0.013 


0.042 


CV Mon 


vdBergh 1 


1-5,7 


0.037 


0.037 


WZ Car 


ASCC 63 


1-5, 7 


0.000 


0.000 



Table 3. Membership probabilities calculated for Combos in- 
consistent with membership. Italic cluster names indicate 'alter- 
native' hosts considered. Constraints by numbers: l=separation 
(i.e. position & cluster radius), 2=7r, 3=RV, 4=/^^/!, 5=fi]jE, 
6=[Fe/H], 7=age. P(A|B) = P(A) ■ P(B|A) , see Sec.[2] 



Cepheid 


Cluster 


Constraints 


P{A\B) 


P(B\A) 


CS Vel 


Ruprecht 79 


1-5,7 


0.148 


0.148 


V1726 Cyg 


Platais 1 


1,3,4,5,7 


0.084 


0.120 


V442 Car 


NGC 3496 


1,4,5 


0.013 


0.035 


TV CMa 


NGC 2345 


1-5,7 


0.000 


0.026 


WZ Sgr 


Sgr OB7 


1,2,3,7 


0.000 


0.012 


SZ Tau 


NGC 1647 


1-5,7 


0.000 


0.010 


UY Per 


Czernik 8 


1,2,4,5,7 


0.000 


0.001 


SU Cas 


FSR 0581 


1,2,4,5,7 


0.000 


0.000 


V1334 Cyg 


ASCC 113 


1,3-7 


0.000 


0.000 


TW Nor 


Nor OB5 


1,2,3,7 


0.000 


0.000 


RU Set 


Dolidze 34 


1 


0.000 


0.000 


SW Vel 


NCC 2660 


all 


0.000 


0.000 


S Vul 


Turner 9 


1-5,7 


0.000 


0.000 


SV Vul 


Turner 9 


1-5,7 


0.000 


0.000 



The remaining 4 bona-fide CCs are not included here, 
since they were found to lie outside the cluster radii defined. 
3 of them are located within two n. V379Cas lies at 4.8 r; 
from its host cluster's center coordinates. 



4- 1-4 Inconclusive Combos and Non-members 

Combos for which our results are inconclusive are listed 
in Tab.Il The Cepheids GSC 03642-02459 (BD-H474225), 
Y Car, and WZ Car have thus far not been considered for 
cluster membership in the literature. 

Two kinds of Combos inconsistent with membership are 
hsted in Tab.[31 Italic cluster names identify Combos for 
which only 'alternative' host clusters were considered, i.e. 
ones whose literature host clusters were not included in the 
cluster sample considered. The remaining entries in this ta- 
ble are Combos considered inconsistent with membership, 
but have been previously considered for membership in the 
literature. 

More details on the Combos listed in Tabs.[2]and[3]are 
provided in the appendix. 



4.2 The Galactic Cepheid PLR revisited 

Let us now employ our bona-fide CC sample in a calibration 
of the Galactic Cepheid PLR. It represents an ideal sample 
to this end, due to the high confidence we can have in cluster 
membership. 

A Cepheid's absolute V-band magnitude. My, is de- 
termined using the true distance modulus of the cluster, 
{Vo — Mv)qi, the Cepheid's mean magnitude, (mv), the ra- 
tio of total-to-selective extinction towards the cluster, -Rv',cij 
and the Cepheid's color excess, E{B — V), as 



Mv = (mv) - (Vo - Mv)ci - Rv,ciE{B - V) , 



(18) 



where quantities refer to the Cepheid, unless subscripted by 

'cr. 

As described in Sec.|3] we compile cluster data from 
the references that yield the largest limiting radii. Obvi- 
ously, this preference (bias) does not necessarily select the 
'best' data available. For instance, line-of-sight dependen- 
cies of reddening are not taken into account by K05 and 
Bll, and the ratio of total-to-selective absorption is fixed at 
the canonical value of 3.1. However, there may be a benefit 
in using the cluster data from these large catalogs, since the 
data are relatively homogeneous. In the following, we refer 
to the set of data compiled for our membership analysis as 
the 'standard' set. 

In an attempt to improve accuracy, we compile data 
from detailed studies of the host clusters. Accurate true dis- 
tance moduli based on ZAMS-fitting for many clusters can 
be found in TIO, for instance. The values given are par- 
tially based on published literature and partially re-analyzed 
using 2MASS and unpublished UBV data (Turner 2012, 
priv. co mmunication) . For N GC 7790, we adopt the detailed 
work by iGupta et all l|2000l ). and for th e two ASCC clus- 
ters, w e use the de-red dened values by iKharchenko et ahl 
(2005aD. For reddening, lAn et al.l |2003) provide a conve- 
nient way to calculate Rv for most lines of sight in our 
sample as Rv = Rv,o + 0.22(5 — V)o, with Rv,o tabulated 
for the cluster and taking into account the intrinsic color 
of the Cepheids. For the two ASCC clusters discovered by 
IKharchenko et all (|_2005a') , we employ Rv = 3.1, since no 
othe r estimate is availab le. For Turner 9 we use Rv = 2.94 
as m [Turner et all (|l997l ). The Cepheid mean V-magnitudes 
and E{B — V) were compiled as described in SecO the 
mean magnitudes were taken typically f rom K S09, and the 
E{B - V) either from iKovtvukh et al.l (|2008| ) or from the 
Fernie database. The values employed are listed in Tab.jlJ 
and we refer to this data set as the 'optimal' one. 

Figure[T6lshows the PLR fits to both the 'standard' (left 
panel) and the 'optimal' (right panel) data sets. In both 
fig ures, four straigh t line s indicate the PL R calibrations 
by Tammann et al.l (|2003l . red dash-dotted) , iTurner et al.l 
(|20ld . green dotted, lower zero-point), as well as our non- 
weighted (dashed cyan), and weighted (solid blue) least- 
squares fits to the data. The large scatter (RMS=0.65) in the 
'standard' data set is striking. Contrastingly, the 'optimal' 
set appears worthy of its name, with an RMS of 0.21 mag. 
We determine the uncertainties on the 'optimal' set by linear 
regression and obtain: 

Mv = - (3.08 ± 0.50) log P - (0.94 ± 0.42) . (19) 
Albeit by chance and with much larger uncertainty. 



Cepheids in Open Clusters 15 



03 




My =-2.96 ±0.69 . log(P) -1.02 ±0.57 
RMS =0.65 

J I I I I I I 



0.5 0.6 0.7 0.8 0.9 1.0 1.1 

log(P [d]) 



-5.0 
-4.5 
-4.0 
-3.5 
-3.0 
-2.5 
-2.0 
-1.5 



1 1 r 


1 1 1 1 

~p ^ 












My =-3.08±0.50 ■ log(P) -0.94±0.41 ~ 


1 1 1 


RMS =0.21 

1 1 1 1 



0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 

log(P \d]) 



Figure 16. Cepheid PLRs fitted to the 'standard' (left panel) and 'optimal' (right panel) data sets. Solid circles indicate previously 
known CCs, green solid circles highlight the three Cepheids in NGC 7790. Two s olid red diamonds in dicate the two new bona-fide 
fundamen tal-mode CCs S Mus an d SX Car. The green dotted line shows the PLR by [Benedict et al.|[2007l . and the red dash-dotted line 
represents iTammann et alj|2003 . Blue solid and cyan dashed lines indicate weighted and non- weighted least-squares fits. An accurate 
PLR calibration critically depends on accurate distance estimates from detailed studies that include line-of-sight- variations of extinction. 



Cepheid 


Cluster 


(mv) 


E{B - V) 


[Vo - Mv) 




logP 




U Sgr 


IC 4725 


6.72 


0.40 ±0.02 


8.81 ±0.10 


3.32 ±0.22 


0.829 


-3.41 ±0.19 


CF Cas 


NGC 7790 


11.15 


0.53 ±0.03 


12.60 ±0.15 


3.33 ±0.27 


0.688 


-3.20 ±0.27 


CE CasB 


NGC 7790 


10.94 


0.50 ±0.05 


12.60 ±0.15 


3.33 ±0.27 


0.711 


-3.34 ±0.34 


CE CasA 


NGC 7790 


11.09 


0.48 ±0.05 


12.60 ±0.15 


3.31 ±0.27 


0.651 


-3.10 ±0.33 


DL Cas 


NGC 129 


8.98 


0.49 ±0.02 


11.11 ±0.02 


3.46 ±0.23 


0.903 


-3.82 ±0.20 


SUCyg 


Turner 9 


6.89 


0.07 ±0.02 


9.33 ±0.20 


2.94 ±0.39 


0.585 


-2.66 ±0.22 


V Cen 


NGC 5662 


6.87 


0.17 ±0.05 


9.28 ±0.05 


3.47 ±0.39 


0.740 


-2.90 ±0.25 


S Nor 


NGC 6087 


6.41 


0.27 ±0.05 


9.78 ±0.03 


3.74 ±0.86 


0.989 


-4.37 ±0.42 


V340 Nor 


NGC 6067 


8.41 


0.34 ±0.01 


11.19±0.15 


3.36 ±0.30 


1.053 


-3.92 ±0.21 


S Mus 


ASCC 69 


6.13 


0.15 ±0.03 


10.00 ±0.20 


3.10 ±0.51 


0.985 


-4.34 ±0.27 


SX Car 


ASCC61 


9.12 


0.33 ±0.03 


11.14 ±0.20 


3.10 ±0.51 


0.687 


-3.04 ±0.33 



Table 4. Parameters adopted for the 'optimal' set used in Eg . 1181 and the right panel of Fig. 1161 True distance moduli of clusters, 
(Vo ~ My), and absorption-relevant parameters were adopted according to the criteria specified in the text. We adopt 0.04 mag as the 
uncertainty on (mv). 



this is nearly identical to the result of Mv = 
- (3.087 ± 0.085 ) log(P ) - (0.914 ± 0.098) obtained by 
iTammann et al.l (|2003l ). Owing to the large uncertainties 
of the fit, our solution is in agreement wit h most other 
zero-p o int-calibrated PLRs such as the ones in lFouaue et al.l 
l|2007l ): [Turner et al.l (12 010'). We do note, however, that the 
solution by [Benedict et al . (200^) lies rather far from our 
calibration in both slope and intercept. 

We note that no correction for binarity was applied in 
this calibration and that all Cepheids were treated as sin- 
gle stars, which many of them are not. The calibration is 
mainly limited by the accuracy of the cluster data, so that 
a correction for binarity is not likely to impact our result 
significantly. 



5 DISCUSSION 

5.1 Constraining Power and Limitations of 
Membership Constraints 

5.1.1 The Prior P( A) 

The form of the prior was motivated by radial density pro- 
files of star clusters, see Sec. l2.1l However, the degree with 



which the distribution of stars in the cluster is known varies 
greatly between clusters, and deviations from circular cluster 
shapes were ignored for internal consistency. Furthermore, 
crowding, great distances (~ kpc) to host cluster candidates, 
differential reddening, sparsity, etc. all conspire to compli- 
cate the definition of cluster radii. It is thus not surprising 
that the prior does not perform extremely well as a member- 
ship constraint if taken at face value. Nevertheless, it does 
help to separate the interesting cases from the majority of 
null matches, since it reduces the question of proximity on 
the projected sky to a single number that contains informa- 
tion on the density of the cluster, since both cluster radii 
are used in its definition. 



5.1.2 The Likelihood P(B\A) 

Intuitively, the greatest set of constraints used provides the 
tightest constraints on membership for any cluster-Cepheid 
combination (Combo). Figure[T7] illustrates this. It shows 
a logarithmic normalized histogram of likelihoods for two 
cases: more than 3 constraints used to calculate P(B|A) 
(red); all constraints used to calculate P(B|A) (blue). How- 
ever, the constraining power of a given set of constraints is 
not merely a function of its size. Here we discuss the mem- 
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Likelihood P(B|A) 

Figure 17. Log-normalized histogram of likelihoods. Red: Com- 
binations with 3 or more parameters used for P(B|A). Blue: Com- 
binations with all membership constraints. The more membership 
constraints are employed, the more separated are high and low 
likelihoods, i.e. the better constraint is membership. 
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Figure 18. Likelihoods computed as a function of proper 
motion. Abscissa: proper motion of Cepheid divided by the 
squared-summed uncertainties, i.e. fi = ^ (/^S)^ 

\l a (/.tj^)^ -(- o" (m^)^- Ordinate; same for Cluster. 
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bership constraining power of the different constraints used 
to calculate likelihoods. 

Proper motions can be very effective at ruling out 
membership if the motion clearly exceeds the uncertainties 
on the measurement. This is the case for a cluster cross- 
matched with a background Cepheid, for instance. However, 
for a significant fraction of Combos, the proper motion vec- 
tor's magnitude was smaller than the uncertainty of the mea- 
surement, thus effectively not constraining membership, see 
FigUHl K the magnitude surpasses the uncertainties by at 
least a factor of 3, proper motion serves as a reliable con- 
straint. For the majority of Combos that fulfill this criterion, 
membership tends to be excluded. 

Distance is a potentially very strong membership con- 
straint, since intuitively, a Cepheid that occupies the same 
space volume as a cluster should be a member. However, 
cluster distances can be subject to large systematic uncer- 
tainties due to parameter degeneracy (distance, age, red- 
dening), model-dependence (rotation, etc.), or previous dis- 
tance estimates to, e.g., the Pleiades. Furthermore, implicit 
assumptions on cluster membership can significantly impact 



the distance determined, especially for relatively young clus- 
ters that harbor few stars around the Main-sequence turn- 
off and few or no red giants. Very detailed studies of open 
clusters, and in particular of the line-of-sight extinction, 
are required for improvement in this domain, as is shown 
in Sec. 14. 21 or demonstrated by the discussion of CVMon's 
membership in van den Bergh 1 (cf. Sec. lAf .0.3|) . 

Radial velocities have the potential to provide very 
tight membership constraints, since the RV dispersion 
within open clusters can be significantly below fkms~^ 
iLovis fc Mavoill200'i1 ). approaching the measurement pre- 
cision on v-y for non-binary Cepheids. However, estimates of 
average cluster RVs are usually based on only a few stars, see 
Fig.[5] fn fact, approximately half of all clusters with 'known' 
average RVs are based on measurements of two stars or less, 
and strong selection effects (e.g. toward late-type stars) can 
severely impact the estimate. Since Cepheids are very bright 
and of late spectral type, it is rather likely that a cluster RV 
is in part based on measurements that include a Cepheid. We 
therefore highlight the need to observe more radial velocities 
of upper main sequence stars in clusters in order to ensure 
the most accurate estimates of average cluster RV. Alterna- 
tively, larger telescopes may observe the much fainter lower 
main sequence . The Gaia-ESO public spectroscopic survey 
ICilmore et all (|2012l ) will soon provide precise RVs for a 
large number of clusters and thus can improve the reliabil- 
ity of a future study similar to the present one. 

The iron abundances compiled here were, arguably, 
of limited use as membership constraints, since data inho- 
mogeneity, the limited number of cluster stars used for de- 
termining the cluster average, and differences in the solar 
reference values used are at the same order of magnitude as 
the range of iron abundances found in the sample consid- 
ered (that lies within the young metal-rich Galactic disk). 
For a few cases, however, the iron abundance did further 
strengthen the interpretation of excluded membership. 

Age as a membership constraint quantifies valid evo- 
lutionary considerations that are established empirically, ft 
is a particularly useful membership constraint, since it is 
readily available for clusters as well as for Cepheids (from 
period-age relations), especially when few other constraints 
are available. However, ages for both kinds of objects are 
subject to model-dependence, and the accuracy of the val- 
ues inferred is difficult to quantify. 

All of the above quantities have their own peculiarities, 
and thus no single one can be named the 'best' membership 
constraint. Instead, the greatest constraining power resides 
in the combination of all the data, as is seen in the hona-fide 
CCs identified in this work, as well as in Fig. ll7l 



5.1.3 Membership Probabilities P(A\B) 

We computed P{A\B) simply as the product of the prior 
and the likelihood, leaving out the normalization term P(B) 
that would in principle be required, see Eq.[T] However, in 
order to make full use of the P{A\B) values computed, P(B) 
should not be neglected. Given that the incompleteness for 
open clusters is quite significant at heliocentric distances 
greater than, say, 1 kpc, we did not currently see this as 
feasible, cf. Fig.[Tl 



Cepheids in Open Clusters 17 




2 4 6 8 10 



d [kpc] 

Figure 19. Distance distribution of known open clusters in log- 
age range [7.0, 8.5] (blue line with high peak at shorter distances) 
and of Cepheids at low Galactic latitudes, 10° (green line). 

Distances of hona-fide CCs from this work indicated by short 
black lines at the top. 



5.2 Incompleteness 

Despite our aim to maximize the number of clusters and 
Cepheids considered, only 13 CCs were identified. This is a 
result of the apparent difference in sample completeness for 
clusters and Cepheids. Fig.[l9] shows the histograms of he- 
liocentric open cluster distances for clusters withing the ap- 
pro priate age range ( logm (iffefyr]) = [7.0, 8.5]), taken from 
the iDias et ali l|2002l) catalog as the blue distribution (left 
peak), and the Cepheid distances compiled (cf. Sec. 13. 2. fl 
green distribution). It is evident that the detection rate of 
open clusters falls off quickly at distances in excess of 1 kpc, 
probably since their identification against the field becomes 
increasingly difficult. Cepheids, on the other hand, are de- 
tectable at much greater distances, thanks to their high lu- 
minosity and characteristic brightness variations. 

However, a cluster's probability of hosting a Cepheid 
is governed by stellar evolution and star formation. For a 
given distribution of stellar masses, only a small fraction will 
become seen as Cepheids during their lifetime (stars with 
masses between ~ 4— IIMq ) . Due to the nature of the IMF, 
these intermediate mass stars constitute a small fraction of 
the total number of cluster stars, and few such intermediate- 
mass stars are present in typical (small) open clusters. Fur- 
thermore, only a fraction (perhaps 10-20per cent) of a suit- 
able star's lifetime during the core Helium burning phase is 
spent on the blue loops, and even less within the instability 
strip. As a result, few CCs are known. The inverse problem 
of finding host clusters around known Cepheids, however, 
can be interesting, e.g. to constrain the survival rate of open 
clusters. 



5.3 Implications for the Distance Scale 

We are aware of the fact that our calibration, although per- 
formed on an 'optimal' set, remains inhomogeneous in the 
true distance moduli employed and the treatment of extinc- 
tion. It is furthermore based on V-band data obtained in 
multiple different passbands with varying post processing 
techniques. These are the main limitations of the data em- 
ployed. Finally, the number of CCs employed in the fit (10) 



is quite small, and the distance to NGC 7790 enters the fit 
three times, since its three CCs are included in the fit. 

We advocate that the 'optimal' sample presented, al- 
though incomplete (missing host clusters), forms an ideal set 
for an open cluster-based PLR calibration, since confidence 
in membership of all clusters was self-consistently evaluated. 
However, to improve upon the calibration of the PLR, a de- 
tailed homogeneous deep photometric study of the host clus- 
ters that includes careful treatment of extinction is required. 
Given the large discrepancies that are found between recent 
Galactic PLR calibrations (and their zero-points), such an 
observational campaign would be very desirable. 



6 CONCLUSION 

Focusing on Cepheids, we have performed an all-sky clus- 
ter membership census. Our analysis considers an up to 
8-dimensional membership-space that includes spatial and 
kinematic information, as well as parent population param- 
eters (age, iron abundance). Although in some ways lim- 
ited by data inhomogeneity and incompleteness, we recover 
most canonical bona-fide cluster Cepheids accessible to us 
and identify three interesting new bona-fide CCs. 

The newly identified CCs are: SXCar in ASCC61, 
SMus in ASCC69 (both fundamental- mode pulsators), and 
V379 Cas in NGC 129 (overtone pulsator). The cluster mem- 
bership of these three Cepheids must have escaped previous 
discovery, since the host clusters are not very well-studied, 
and the Cepheids are located outside the cluster cores, cf. 
Sec.O 

Based on the available data, we find our results to be 
inconsistent with membership for the following combina- 
tions previously considered in the literature: V442 Car in 
NGC 3496, TVCMa in NGC 2345, V1726Cyg in Plataisl, 
UYPer in CzernikS, WZSgr in SgrOB7, SZTau in 
NGC 1647, and CS Vel in Ruprecht 79. 

Since we can rank candidates through membership 
probabilities, we consider our bona-fide CC sample to be 
ideal for calibrating the Galactic Cepheid PLR. In Sec. 14. 21 
we perform this calibration and determine My = —3.08 ± 
0.50 log P — 0.94 ± 0.42 to be our best sol ution, which is 
in cha nce agreement with the calibration bv lSandage et af] 
However, this calibration remains limited in preci- 
sion, due to lack of homogeneity and rather uncertain clus- 
ter data. We therefore highlight the need for observational 
campaigns dedicated to the host clusters of our bona-fide 
CC sample. 

The limitations that our work suffers due to inhomoge- 
neous and incomplete data will be significantly reduced in 
the near future, thanks to the Gaia space mission. Specifi- 
cally, Gaia will improve our study in the following ways: 

• Thousands of ne w Cepheids (|Ever fc Cuvperd I2OO0I : 
IWindmark et al.ll201ll ) will be discovered. 

• Accurate absolute trigonometric parallaxes of Cepheids 
will become av ailable up to dista nces of 6-12 kpc, depending 
on extinction (|Ever et al.ll20lj ). thereby enabling a direct 
calibration of the Galactic period-luminosity-relationship 
similar to the on e perf ormed in the seminal paper by 
iFeast fc Catchpold l| 19971 ). This will remove our partial de- 
pendence on existing PLR calibrations when determining 
cluster membership. Accurate parallaxes to longer period 
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Cepheids will also be obtained, thereby significantly improv- 
ing the distance estimates to extragalactic Cepheids. 

• The accuracy of proper motions will be improved by 
orders of m agnitude, moving from masyr"^ to tens of 
fiasyr~^, see iLindegrenI (|2010l ) and the Gala Science Per- 
formance websit^l 

• Homogeneous radial velocities (via the RVS instru- 
ment) and metallicity estimates (via the spectrophotometric 
instruments) will be available as membership constraints for 
a great number of Cepheids. 

• Homogeneous metallicit y estimates can be obtained 
through spectro-photometry l|Liu et al.l 20121) . 

• Thousands of new open clusters (|ESA| I2OO0I ) will be 



discovered, forming a more or less complete census of open 
clusters to distances up to 5 kpc. As was shown in Fig. 1191 the 
distribution of Cepheids is still increasing at these distances. 
Hence, one may expect to find many more CCs in these parts 
of the Galaxy. 

• Down to magnitude 20, all-sky homogeneous multi- 
epoch photometry and colors will be obtained that will in- 
clude all the bona- fide clusters mentioned in this work. 

• Known clusters will be mapped in unprecedented detail, 
and intra-cluster dynamics will be accessible to determine 
membership. 

Gala's data homogeneity will significantly improve error 
budgets, since no offsets in instrumental zero-points (e.g. 
in RV) will have to be taken into account. The constraining 
power in terms of membership will thus be augmented con- 
siderably. Correlations between parameters, e.g. proper mo- 
tion, parallax, and R V, can be determine d self-consistently 
and accounted for (cf. Ivan LeeuwenlboOTl V Such factors will 
make the Gala era particularly exciting for work such as 
this. 
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APPENDIX A: DETAILS ON INDIVIDUAL 
CLUSTER CEPHEID COMBINATIONS 

Al Inconclusive CC Candidates 

Al.0.1 CG Cas - Berkeley 58 CG Gas lies at a sepa- 
ration of 5.7' (outside the core) from Berkeley 58 's center, at 
roughly half the limiting radius. No average cluster RV was 
found. While the Gepheid's PLR-based parallax is close to 
that of the cluster and reddening is in agreement, age tends 
to contradict membership. Furthermore, the differing proper 
motion (p* ) does not suggest a common point of origin with 
the cluster and therefore strongly disfavors membership. In 
light of this, we cannot consider GG Gas a bona-fide cluster 
Cepheid. A detailed investigation of cluster RV and proper 
motion would be helpful to clarify the membership status. 

Al.0.2 GHCar and Trumpler 18 Membership of 
the binary overtone Cepheid GH Car jSzabados et al.l 
20121) in cluster Trump ler 18 was first proposed by 
Vazquez & Feinstein (199C|) and then called into question by 
Baumgardt ct al. (2000,) who recommended radial velocity 
follow-up to draw a firmer conclusion. We compute a high 
likelihood of 57 per cent based on proper motion, age, and 
RV. No parallax was calculated, since GH Car is an over- 
tone pulsator. However, the Gepheid's distance listed in the 
Fernie database (2.2 kpc) is significantly larger than Trum- 
pler 18's (1.4 kpc). Furthermore, the Cepheid is reddened by 
0.1 mag more than the cluster, which is consistent with a 
greater distance to the Cepheid. 

Trumpler 18 is sparsely populated (17 members consid- 



ered in K05), and there seems to be some uncertainty regard- 
ing its age, as seen by the great age difference in K05 and 
D02 (the latter would be inconsistent with the Gepheid's 
pulsational age). The average cluster RV is based on the 
measurements of a single star, and does not agree with that 
of the Cepheid. Cepheid and cluster are co-moving in proper 
motion, however, with only a small difference in fis- In sum- 
mary, there is considerable doubt regarding GH Car's mem- 
bership in Trumpler 18. A more detailed study of the clus- 
ter's parameters, including its average RV would help to 
clarify this situation. 



Al.0.3 CVMon and van den Bergh 1 CVMon lies 
right in the center of cluster van den Bergh 1 and was 
previously considered to be a member of this cluster, cf 
Turnerj3t_al. (1998). The constraints employed in our cal- 
culation are parallax, proper motion, radial velocity, and 
age. From these, we compute a likelihood of merely 4 per 
cent. The cluster data from Bll is largely inconsistent with 
the data for the Cepheid, see e.g. parallax (also redden- 
ing) and age. The average RV given in D02 is identical 
to that of CVMon and no information on the number of 
stars e mployed for its determin ation is given in the original 
source (|Rastorguev et alll 19991 ): it should thus be discarded 
as membership constraint. The same goes for the cluster's 
proper motion which w as determined on CV Mon alone by 
iBaumgardt et al.l (|2000l ) using the original Hipparcos cata- 
log (hence an offset exists to the proper motion of CV Mon 
in our data compilation). 

[Turner et al.l ((1998') addresses a number of issues en- 
countered when studying this cluster, such as the possible 
existence of oth er clusters very nearby van den Bergh 1. 
lAn et al.l |20o3) also noted difficulties involved with this 
particular cluster, and there may be unreddcned field stars 
present that contaminate the Main Sequence (Pres t onlll963 ) 
impacting d istance and age estim ates. The cluster's distance 
estimate by [Turner et al.l l| 19981 ) is in close agr e ement with 
the Gepheid's IRSB-distance bv iFouaue et al.l (|2007l ). and 
their cl uster reddening very c lose to the Gepheid's E(B- 
V) from iKovtvukh et al.l (|2008h . Furthermore their average 
cluster RV determined on upper main sequence stars also in 
good agreement with CVMon's v^. 

To summarize, two of the four membership constraints 
compiled are not valid for this particular case, and the de- 
termination of distance & age is difficult. The issues en- 
countered here are due to a lack of detail on cluster param- 
eters in the catalogs, especially for a complex case as this, 
and demonstrate that detailed studies focusing on individual 
clusters continue to be of great importance. 



Al.0.4 BE Sgr and Collinder 394 BB Sgr lies at 21' 
separation from Collind er 394's center, and membership was 
first studied in detail bv lTurner fc Pedreros! (Il985l ) after hav- 
ing been originally suggested by Tsarevskv et al.ll|l966l ) . Un- 
fortunately, the data available as membership constraints 
neither really supports, nor excludes, it: The mismatch in 
RV, although large, is absorbed by the relatively large un- 
certainties, and parallax, reddening, and age agree. The only 
real discrepancy is seen in ^* , although also this discrepancy 
barely exceeds the mutual uncertainty. 
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Al.0.5 Y Sgr and IC 4725 (M25) The parallaxes em- 
ployed (Cepheid ttcop = 2.13±0.29 from lBenedict et al.l2007l 
and Cluster ttci = 1.61±0.32) in the calculation nearly agree 
within their respective error budgets. However, the lower 
reddening of the Cepheid may make Y Sgr seem to lie in the 
foreground of M 25. 

Using the well-established CC U Sgr as a point of ref- 
erence, we remark that proper motion, age, and metallicity 
are in excellent agreement between both Cepheids. v-y differs 
slightly between the two, which could be explained by the 
known binarity of USgr and YSgr. The only discrepancy is 
in distance, which might be explained by uncertainties in ex- 
tinction, since M25 lies in the Orion arm (XHIP). Although 
we calculate a very low prior for this Combo, Y Sgr's distance 
from cluster center is merely 4.2 pc (projected to M25's dis- 
tance). From these considerations, it appears that YSgr is a 
could be a member candidate, pending confirmation through 
follow-up of extinction. 

Al.0.6 BD+47 4225 (GSC 03642-02459) and 
Aveni-Hunter 1 The star was cla s sified as a Cepheid 
based on HAT data by iBakos et all (|2002l ') with a light 
curve that suggests fundamental-mode pulsation. It lies 
barely outside Vc of cluster Aveni-Hunter 1. Unfortunately, 
proper motion and age are the only available membership 
constraints available in the literature compiled. Cluster and 
Cepheid appear to be co-moving in proper motion. Further- 
more, the pulsational age of the Cepheid is spot-on with the 
cluster. However, a very rough distance estimate using the 
V-band magnitude (10.47) found in the Guide Star Catalog 
V. 2.3.2 (|Lasker et al.|[200^ yields a distance of 2.5 kpc for 
the Cepheid, which is 5 times the cluster distance. Since 
both objects are thus far not very well-studied, we highlight 
the need for follow-up of both cluster and Cepheid. 

Al.0.7 YCar and ASCC 60 YCar is a double-mode 
Cepheid in a triple system with a B9. 0V companion 
(|Evanj 1 19921 ') on a known orbit (see the ISzabadosI |2003| 

binary Cepheids databasf^. The Cepheid lies well inside 
ASCC 60's projected core. Since the RV of the cluster in 
K05 was measured on a single star and is identical to Y Car's 
v-f, we cannot consider this a valid membership constraint. 
The age estimates lie just outside each other's error bud- 
gets, and reddening agrees. However, the Cepheid's clearly 
discernible proper motion points to a shorter distance, how- 
ever, and is in excellent agreement with the proper motion 
of the cluster. A detailed photometric and radial velocity 
study of ASCC 60 is required to conclude on the possible 
membership of YCar. 

Al.0.8 WZ Car and ASCC 63 Due to the large sepa- 
ration from ASCC 63's core (32') and starkly differing line- 
of-sight velocities (5frad ~ 120kms~^), membership would 
be excluded for this combination, if the probabilities were 
to be taken at face value. 

Looking at the individual constraints, however, we find 
that age, reddening, parallax (IRSB), and proper motion 
strongly suggest membership. Suspiciously, the cluster RV 
is based on only 2 stars, and may therefore not be reliable, 

1* http : //www . konkoly . hu/CEP/intr o . html 



or point towards an ejection event. In any case, this com- 
bination remains interesting and follow-up of the cluster is 
warranted. 

Al.0.9 VZ CMa and Ruprecht 18 This combination 
yields a likelihood of 83 per cent, since metallicity and 
age are in excellent agreement between the cluster and 
the Cepheid. Proper motion is better discernible for the 
Cepheid than for the cluster, it seems, and reddening for 
the Cepheid is less strong than for the cluster. The cluster 
(1.7 kpc) is a fair bit farther than the Cepheid (1.2 kpc, from 
the Fernie database). VZCMa may lie in the foreground of 
Ruprecht 18. However, it would be beneficial to redetermine 
the Cepheid's distance, e.g. with an IRSB technique, and 
obtain an average RV for the cluster to better constrain 
membership for this Combo. 

A2 Combos Inconsistent virith Membership 

As is clear from the results of our analysis, most cross- 
matched Combos are inconsistent with membership. We feel 
it is not necessary to discuss all of these here, and refer the 
reader to the electronic table. However, a few Combos merit 
a short description, for instance since they have been con- 
sidered for membership previously, or lie close to the center 
of a cluster. These Combos are discussed in the following 
subsections. 



A2.1 Previously Discussed in Literature 

A2.1.1 SZTau and NGC 1647 The membership of 
SZTau in the halo of NGC 1647 was first consid ered by 
lEfremovl (| 19641 ') and later studied in more detail bv iTurnCT 
(|1992| ) who concluded that membership is likely. Almost all 
membership constraints could be applied for this combina- 
tion. It is conspicuous that parallax and age agree well for 
cluster and Cepheid (K05), although these are subject to 
model-dependence and the cluster's age listed in D02 would 
not agree with SZ Tau's pulsational age. The radial velocities 
agree within their prescribed uncertainties (5wrad = 6kms~^. 
Based on proper motion and the very large separation of 
more than 2°, however, members hip can be effectiv e ly ex - 
cluded, as was noted previously bv lBaumgardt et al.l (|2000l ) 
using (the old reduction of) Hipparcos proper motions. 

A2.1.2 V1726 Cyg - Platais 1 V1726 Cyg's member- 
ship results for Platais 1 are based on separation, proper 
motion radial velocity, and age. Th e star w as fi rst considered 
for clu ster membership bv |Platai3 |l97^ and [Turner et af] 
(Il994h . In the available constraints, the only real match be- 
tween the cluster's and the Cepheid's quantities is in ra- 
dial velocities. We did not calculate a PLR-based paral- 
lax, since V1726 Cyg is an overtone pulsator. However, the 
Fernie database places V1726Cyg much closer (2 kpc) than 
Platais 1 (3.5 kpc, K12), which is consistent with the signifi- 
cantly lower (0.54 mag) reddening of the Cepheid. fis nearly 
vanishes for the cluster (K12), while V1726Cyg has clearly 
detected proper motion in this direction. Finally, the clus- 
ter is significantly older than the Cepheid. Taken together, 
these constraints yield a very low likelihood that indicates 
non-membership. 
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A2.1.3 CSVel and Ruprecht 79 Membership 
of C S Vel in Ruprecht 79 wa s th oroughly discussed 
by [ Harris fc van den Berghl l| 19761 ) who credited 
iTsarevskv et alj (Il966l ) with first suggesting tliis par- 
ticular combi nation. It has since been studied multiple 
times, e.g. by IWalkeJ (|l987l l and TIO. Due to the sparse 
nature of the cluster, i ts reality as such was doubted by 
iMcSwain fc Gied (12005'), who conclude that Ruprecht 79 
rather be a hole in the dust of the Sagittarius-Carina spiral 
arm than a physical open cluster. 

The data compiled in this work would place the Cepheid 
about 1 kpc farther from the Sun than Ruprecht 79 (with 
E(B-V) by 0.1 mag higher for the Cepheid than the cluster), 
while the average cluster RV (K05) is not determined clearly 
enough to provide a strong constraint. We furthermore note 
that proper motion for the cluster and the Cepheid strongly 
contradict membership. Therefore, even if Ruprecht 79 is, in- 
deed, an open cluster, the information available would be in- 
consistent with membership, yielding a likelihood of merely 
5 per cent. 

A2.1.4 V442 Car has pr eviously been co n sidere d for 
membership in NGC 3496 by iBalona fc LanevI (119951 ) who 
concluded it to be a background star, based on age and red- 
dening. From proper motion and separation, we come to the 
same conclusion. Furthermore, a rough distance estimate for 
a 14"^ magnitude 5.5 d Cepheid excludes membership in a 
cluster located approx. 1 kpc from the Sun. 

A2.1. 5 UYPer and CzernikS or King 4 iTurneij 
(| 19771 ) suggested that UYPer could be a member of Cz- 
ernikS or King 4. [Turner et al. I dioTo") again mentions the 
latter combination. Our results, however, are inconsistent 
with membership in either cluster, based on the constraints 
parallax, proper motion, and age. 

The 'likelier' of the two Combos is CzernikS, for which 
parallax and age are in relatively good agreement; the 
Cepheid is slightly farther away and has larger reddening. 
However, kinematically, the cluster's vanishing proper mo- 
tion (K12) would be inconsistent with membership of the 
rather fast moving Cepheid (/i* — —6.15 ± 2.Smasyr~^, 
fis = 12.S9±2.9masyr-\ PPMXL). 

A2.1.6 TVCMa and NGC 2345 The membership 
constraints parallax, proper motion (vanishes), and age, 
agree within their respective uncertainties. Radial velocity 
differs by ~ 20kms~^ between cluster and Cepheid. In line 
with the prior, we take this to indicate non-membership for 
this Combo. 



A2.2 Foreground Cepheids 

A2.2.1 OPCas and AU Cas in Cas OB7 Although 
the parallaxes compiled for OP Cas and Cas 0B7 agree 
within their uncertainties, the large difference (~ 1 mag) 
in reddening is indicative of non-membership. No Cepheid 
RV is available, and proper motion vanishes. Although the 
age agrees well, we cannot have confidence that the cluster 
parameters are correct, e.g. if reddening was overestimated. 
The same is true for AU Cas. 



A2.2.2 COAur and Koposov 12 All available in- 
formation points towards a foreground Cepheid: the liter- 
ature distance to CO Aur is significantly shorter than to 
Koposov 12, and proper motion is non-vanishing for the 
Cepheid, while it vanishes for the cluster. 

A2.3 Background Cephetds 

There are a few Combos worth mentioning here that yield 
high likelihoods (> 75 per cent) and low priors (~ per 
cent). These ones are WXLac and IC 1442, CM Cas and 
IC 1S05, as well as XVul and Turner 9. The cluster and 
Cepheid parallaxes for each of these Combos agree within 
the stated uncertainties. However, for all of these Combos, 
the Cepheid reddening is significantly (at least twice, or 
0.65 mag) larger than for the cluster. Larger reddening is 
indicative of greater distance. We therefore dismiss these 
candidates, taking into account also the remaining member- 
ship constraints. 

A2.3.1 V410 Vul and NGC 6823 Located outside 
f'lim, proper motion vanishes for the Cepheid, but not for 
the cluster. The faintness of Cepheid (lS*'^mag) hints at a 
location in the background of NGC6S23. 

A2.3.2 NSV1995T and NGC 5281 The small am- 
pli tude Cepheid (given 2"^^ rank as a Cepheid candidate 
by ICaldwell et al.l Il99ll ) lies outside rc in the vicinity of 
NGC52S1. Although proper motion yields a high likelihood 
(0.S9) and the separation is rather small in absolute terms 
(6.9'), the Cepheid is likely a background object. The mean 
I-band magnitude of 13. IS appears too faint for NGC 52Sl's 
distance of just above 1 kpc. We conclude that membership 
is unlikely, but suggest follow-up of the Cepheid to better 
understand its nature. 

A2.3.3 NSV 19262 and Ruprecht 10 1 The star was 
deemed a good Cepheid candidate by Cald well et al.l (|l99ll ) 
and lies at a separation of 4.4' > rc from Ruprecht lOl's 
center. The faint I-band magnitude of approx. 13.15 together 
with a period of roughly 15 d (estimated from Fig. 7ff in the 
original report) expose the star's nature as a background 
object to the 2.3 kpc-distant cluster. 

A2.3.4 NSVS 2089061 and NGC 1496 Cepheid and 
cluster appear to be co-moving and have matching ages, with 
the Cepheid located in the cluster's limiting radius. More 
information on the Cepheid is required to properly assess 
membership. However, from a rough distance estimate it 
would seem likely that the Cepheid is located at around 3 
to 5 kpc, whereas the cluster is much closer (1.4 kpc). 



